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Abstract 



The Ceo molecule and its fragments are studied using Configuration Interac- 
tion (CI) and Variational Monte Carlo (VMC) techniques, within the Hubbard 
model. Using benzene as a test case, we compare the results of the approxi- 
mate calculations with exact calculations. The fragments of Ceo studied are 
pyracylene, fiuoranthene and corannulene. The energies, bond orders, spin- 
spin and charge-correlation functions of these systems are obtained for various 
values of the Hubbard parameter, U. The analysis of bond orders and corre- 
lation functions of these individual molecules allow us to visualise pyracylene 
as a naphthalene unit with two ethylenic moieties and fiuoranthene as weakly 
bridged benzene and naphthalene units. Corannulene is the largest fragment of 
Ceo that we have studied. The hexagon-hexagon(h-h) bond orders are slightly 
larger than those of the hexagon-pentagon bonds (h-p), a feature also found in 
other fragments. We also find bonds between two coordinated carbon sites to 
be stronger than bonds involving three coordinated carbon sites. In Ceo, the 
h-h bonds are stronger than in corannulene and the h-p bonds weaker than 
in corannulene for all correlation strengths. Introducing bond alternation in 
the buckyball enhances this difference. 
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1 Introduction 



The discovery of Cgo, a new allotrope of carbon, by Kroto et al[]I[ and its 
bulk synthesis by Kratschmer et al[§ have been amongst the most exciting 
recent developments in chemistry and physics. Cqq and its compounds show 
exotic physical and chemical properties. Alkali doped compounds of Ceo are 
superconducting, with -R&i.sCsi.sCeo having a superconducting transition tem- 
perature of 32 K|Q|. The compound Cqq — TDAE {T D AE=tetTakis dimethyl 
amino ethylene) is an organic ferromagnet with a fairly large Curie tempera- 
ture of 16.5K0, ||. 

To develop models for diverse phenomena such as these, it is essential to 
have a reliable understanding of the electronic structure of the basic building 
blocks of these systems. The carbon atoms in fullerenes are in a nearly sp"^ 
hybrid configuration with a half-filled p7r-like orbital directed radially and in- 
volved in extended conjugation. Electronic structure studies on conjugated or- 
ganic systems have brought out clearly the importance of electron correlations 
to account for many of the observed spectroscopic and related properties[^]. 
Explicitly including electron correlations is a formidable challenge even at the 
level of a single buckyball. It is reasonable to hope that a proper study of 
a single buckyball and its fragments will shed light on the esoteric electronic 
properties of these systems in the bulk. 

It would be interesting to try and understand the properties of the full 
buckyball in terms of its fragments. Firstly, for small enough fragments, it is 
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possible to perform exact calculations that allow us to systematically check the 
approximate methods. The electronic structures of the fragments are expected 
to be useful in interpreting the electronic structure and geometry of Cqq. In 
this article, we study the ground state properties such as bond orders, charge- 
charge and spin-spin correlation functions of some fragments of Ceo using 
configuration interaction (CI) techniques and the Variational Monte Carlo 
(VMC) method. We report the results of our VMC calculations on the full 
buckyball within the Hubbard model, with and without bond-alternation. 

The simplest description of the electronic structure of Cgo is found within 
the framework of the Huckel picture, wherein a non-zero transfer integral is in- 
troduced between the p.,^-\ike orbitals on nearest-neighbour atoms(sites). The 
resulting molecular orbitals (MOs) have a high degree of degeneracy due to 
the icosahedral symmetry of the molecule. Neutral Ceo has a five-fold degen- 
erate set of highest occcupied molecular orbitals (HOMOs) and a three-fold 
degenerate set of lowest unoccupied molecular orbitals (LUMOs). The high 
degeneracy of the MOs has been the focal point of the discussion of the elec- 
tronic structure of these systems. For instance, in a first approximation, the 
degeneracy of the LUMOs is considered essential to rationalize the existence 
of ferromagnetic exchange between two monoanions of CeolQ. 

The minimal interacting description for Ceo is given by the single band 
Hubbard model which incorporates only an on-site electron-electron repulsion 
U besides the one-electron part given by the Huckel picture. Drastic truncation 
of the electron repulsion to an on-site interaction is justifiable in a metal where 
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the mobile conduction electrons lead to short Debye-Huckel screening lengths. 
In molecules it is strictly necessary to include extended range interactions, at 
least at the level of the PPP model, because of inefficient screening. However, 
obtaining reliable results even for such models is difficult and we study the 
simpler Hubbard model as a first step in understanding crucial qualitative 
effects of correlations in molecular systems like Cqo and its fragments. Exact 
diagonalization techniques are useful for small systems (with less than about 14 
sites at present) but are impractical for large systems because of the exploding 
dimension of the basis with increasing system size]^]. Therefore, approximate 
approaches to the problem are inescapable for larger systems. However, exact 
diagonalization provides a strong check on any novel approximation scheme. 

Amongst approximate techniques, perturbation theory gives fairly reliable 
results at limiting values of the model parameters but breaks down for inter- 
mediate but relevant values of the parameters]^]. Another equally widely used 
approximate approach has been the variational approach with the Gutzwiller 
wavefunction (GWF) as the trial function of choice for Hubbard models[|T^]. 
In this approach, a variational parameter g {0 < g < 1) is introduced and the 
weight of a configuration in the site representation of the Hartree-Fock ground 
state is modified by a factor g^, where D is the number of double occupancies 
in that configuration. Even though a single variational parameter {g) is to 
be determined, the method itself is not easy to implement. Evaluation of the 
expectation value of the Hamiltonian for the GWF is made difficult by the 
very large number of configurations appearing in the wavefunction in the site 
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representation. Often, uncontrolled approximations are resorted to for eval- 



uating the expectation value |]I0|, [Tl|. Such approximations can be avoided if 
we carry out a Monte Carlo integration for the expectation value, using the 
VMC method. The VMC technique uses an importance sampling of the vari- 
ous Slater Determinants (in the atomic orbital (AO) basis) that make up the 
Gutzwiller trial function to estimate the expectation value of the Hamiltonian. 

We also employ a novel configuration interaction (CI) technique that con- 
siders all many-body states derived from multiple particle-hole excitations 
from the Hartree-Fock ground state, whose one particle energies lie below a 
chosen cutoff. The technique gives reliable ground and excited states of small 
fragments of Cqq for comparison with Monte Carlo studies. We carry out VMC 
studies on the full buckyball using the optimized Monte Carlo parameters ar- 
rived at from such comparisons. 

This article is organized as follows. In the next section we give a descrip- 
tion of the computational schemes we have employed, namely the approximate 
CI scheme and the VMC method employing Gutzwiller trial function. We pro- 
vide detailed comparisons of these techniques with exact results for benzene. 
Section 3 deals with our results and discussions on fragments of Ceo- In sec- 
tion 4 we discuss our VMC results on Cgo and compare these with results on 
fragments. The last section summarizes the main results of the paper. 



6 



2 Computational Schemes 

The computational schemes that are described in this section are apphed to 



the Hubbard model |112| on the buckyball and its fragments. The Hubbard 



Hamiltonian may be written as 

Hnub = Hq + Hi (1) 
^0= + aj^flia) (2) 

<ij> o- 

Hi = Uj2 ^i^^ii (3) 

i 

where Hq is the non-interacting part of the Hamiltonian and Hi represents 
the on-site Coulomb repulsion. The operator a-^ (aia) creates (annihilates) an 
electron with spin a in the orbital at the i^^ site, the corresponding 

number operators and the summation < ij > is over bonded atom pairs. 
The transfer integral in all the computations is fixed at 1.0 and energies are 
reported in units of t unless specified otherwise. 

2.1 VB Method in the MO Basis 

In carrying out approximate electronic structure calculations of large systems, 
within a restricted CI scheme, it is important to find a criterion for truncating 
the basis. In general, it is difficult to find a truncation scheme to restrict the 
size of the Hilbert space, that yields accurate spectral gaps. The crux of the 
problem is that increasing the dimensionality of the basis set in the restricted 
CI scheme undoubtedly improves the ground state energy in accordance with 
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the variational theorem, but does not guarantee a concommitant improvement 
in the excited state energies. The scheme described here is intuitive and does 
not have the drawback of the uncontrolled nature of the multi-reference CI 
(MRCI) schemes in which the MOs chosen in the reference determinants are 
arbitrary||T3|. 

In our restricted CI scheme, we employ the MOs to construct the VB 
functionsiOl. The Hubbard Hamiltonian in the MO basis can be written as 



H = Y. ep^JaV + 1/2 E Wpgrs{EpgErs - SqrEps) (4) 
pa pqrs 

Epq = Y.blbq^ (5) 

a 

where (bpa) corresponds to creation (annihilation) of an electron in the p^^ 
molecular orbital with spin a, is the energy of the p^^ MO. Wpqrs is defined 
by 

^^pqrs U ^ ^ CipCiqCiyCig (6) 
i 

where Cjp is the coefficient of the i^^ AO in an AO expansion of the p^^ MO 
and all the summations run over the entire respective bases. The Hubbard 
Hamiltonian written in the MO basis has the property that the one-electron 
part is diagonal, while the interaction part is off-diagonal. We use the Rumer- 



Pauling||T5[ rules to construct the VB diagrams from the MOs. However, if 
an MO state is degenerate, the associated MO functions may not be orthog- 
onal. To avoid the 'nonorthogonality catastrophe', we construct the Lowdin 
orthogonalized M0s[|l6| whose coefficients are given by 

= S-i/^c (7) 



where S is the overlap matrix in the MO basis and c and are the co- 
efficient matrices. We then form the matrix representation of the Hamilto- 
nian in the MO-VB basis using rules similar to those developed for AO-VB 
calculations |T7|, |T^. The restricted MO-VB basis includes all configurations 



whose energy with respect to the energy of the ground state configuration is 
below a given threshold, Ec- We have observed that when Ec equals the spread 
in the one-particle spectrum (bandwidth of the extended systems), the exci- 
tation gaps smoothly follow the exact gaps in small systems. The limitation 
of this approach arises from the fact that as the system size increases, the 
Hamiltonian matrix becomes less sparse and this places a restriction on the 
largest system size accessible by this technique. If we are only interested in 
the ground state properties of the system, keeping the largest possible cutoff 
would be appropriate in accordance with the variational principle. We usually 
keep Ec ~ bandwidth for ground state properties of a single molecule. 

2.2 Variational Monte Carlo with the Gutzwiller Trial 
Function 



The VMC method employs a trial wavefunction, usually of the Jastrow form|19 



20 1 and provides a Monte Carlo estimate of the expectation value of the Hamil- 
tonian and other operators in this state for different values of the variational 
parameters. Here, we focus on the Gutzwiller trial wavefunction which is the 
simplest Jastrow function and is appropriate for Hubbard models. In this sec- 
tion, we briefly describe our implementation of the VMC method, which is 
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similar to that of Yokoyama and Shiba|^ . 
The Gutzwiller wavefunction is given by 

I^G (8) 

^ = EA (9) 

i 

Di = riiinii (10) 

where g is the sole variational parameter. Here |$ > is the ground state , 
given by, 

1$ >= 1$^ > |$| > (11) 

i<^'.>=n^j.io> (12) 

p 

where the product is over the MO's occupied with electrons of spin cr. The 
state 1$ > can be written in terms of the states |-R >, where \R >= \R) > 
\R^ > represents a given configuration of up and down spins in real space 
(equivalently in AO or site representation), as 

1^ >=j:T.det{C^{Ri))det{C^{RmR^ > \R^ > ■ (13) 

Here det{C\RD) is an x A^^ determinant whose elements in the j*'' row 
correspond to the coefficients of the j*^ occupied MO (for this spin) at the sites 
occupied in the configuration >. The operator D picks out the number 
of doubly occupied sites d{R) in every configuration \R > appearing in |$ > 
(eqn. 13). 
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In the VMC method, we seek to evaluate the energy expectation value of 
the Hamiltonian in the state >, 

- < iik> 

To do so, we rewrite E{g) as 

(15) 

where Pg{R) is the probabilty of occurance of the configuration \R > given by 

R' 

From eqn. (13) for |$ >, it follows that 

< Rnl^G >= g''^''-Uet{C\R\))det{C\Ri)). (17) 

A direct importance sampling of the configurations >} with weights 
given by eqn. (16) is not possible since there is an explosion in the number of 
configurations for large systems and the denominator in eqn. (16) cannot be 
evaluated explicitly. Therefore, wc resort to an indirect importance sampling 
technique which is achieved by a Markov process. The states of the Markov 
chain are the configurations and the one-step transition probability between 
configurations is such as to yield Pg{R) (eqn. (16)) in the crgodic limit. In 
the computational procedure, we generate the configuration \Rn+i > from 
\Rn > by hopping an electron of randomly chosen spin, from a randomly 
chosen occupied site, to a randomly chosen vacant site. The new configuration 
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\Rn+i > is accepted with probability 



|2 



^""^^^ = I < Rnl^a > ? ^''^ 
If the new configuration is not accepted, the configuration |i?„ > is itself 
treated as the new configuration. The ratio Pn.n+i is given by 



,\det{C^{Ri^,))\'\det{C^{RU,))\ 



2 



\det{a{RlW \det{C^{RkW 
where 6d (= d{Rn+i) — d{Rn)) denotes the change in the number of doubly 
occupied sites brought about by the electron hop. We use the updating algo- 
rithm of Ceperley et al|22| to compute the ratio of determinants appearing in 
eqn. (19). This involves updating the inverses of the matrices C"{R") which 
can be carried out in A^^ arithmetic operations^ . 



The energy expectation value is computed as a sum of the transfer part 
and the interaction part. It is trivial to compute the expectation value of 
the interaction part of the Hamiltonian, because it only involves counting the 
number of doubly occupied sites in a given state. The effect of the hopping 
operator Af on a state |-R^ >, that hops an electron from site i is 

Ar|K>=5:iK'^-^> (20) 

i 

where \Rn^^^ > denotes a configuration generated from \R^ > by hopping 
an electron of spin cr from site i to site j. The non-interacting part of the 
Hamiltonian is the sum of over all sites i and spin a. Thus, we can write 
the contribution of the state \Rn > to the total energy expectation value for 
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the Gutzwiller wavefunction as 





\H\Rn > 


< ^gI 


\Rn > 



-'EEE <^ (21) 

The calculation of other expectation values is similar to the calculation of 
energy. In particular, the expectation value of any operator A can be estimated 
by estimating < ^(3|74|i? > / < ^g\R >■ As before, the averaging is carried 
out after allowing sufficient number of moves per electron (MPE) to reach the 
ergodic limit. 

The standard deviation in the estimated quantities can be taken as an er- 
ror estimate in any Monte Carlo procedure. In the estimates of energies and 
related quantities, this requires computing the expectation value of where 
H is the Hamiltonian. The number of resultant Slater Determinants when the 
Hamiltonian acts on a single Slater Determinant is of 0{p), where p is the 
number of transfer terms in the Hamiltonian. In Cqo, P = 360 as there are 
90 nearest neighbour bonds. Therefore, in the estimation of the expectation 
value of H^, the number of Slater determinants resulting from the operation 
of H'^ on a single Slater determinant is of 0{p'^). Thus, computationally, it is 
very expensive to estimate the expectation value of and hence the standard 
deviation. However, small systems are amenable to exact Gutzwiller calcula- 
tions and we have compared VMC estimates with these results to arrive at a 
satisfactory number of Monte Carlo steps. The actual scaling of Monte Carlo 
errors with size of the physical system is not known. In large systems, for some 
model parameters, we have checked that the Monte Carlo averages obtained 
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with 10^ MPE, after allowing about 2000 MPE for reaching ergodicity, do not 
change within the significant figures quoted in the tables even upon doubling 
the sample size. 

2.3 Comparison of Approximate and Exact Results for 
Benzene 

We first compare our results from the VMC method with exact calculations 
in the case of benzene. The calculations on benzene include computation 
of expectation values obtained from the variation of g in the full Gutzwiller 
wavefunction and also a Monte Carlo estimate of the same. We refer to the 
former as the true Gutzwiller value and to the latter as the VMC value. In 
table(l), we compare the energy from exact-diagonalization with the true 
Gutzwiller energy and the VMC as a function of the Hubbard parameter, U. 
We observe that even upon increasing U, both the true expectation value in 
the Gutzwiller wavefunction and the MC estimate continue to compare well 
with the exact energy. 

The ground state energy alone is not sufficient for a physical picture of the 
system and it is essential to have reliable correlation functions. Therefore, we 
have studied the spin-spin and charge-charge correlation functions using all 
the techniques mentioned above. In table(2), we present the spin-spin and 
charge-charge correlation functions obtained from the exact calculation, the 
true Gutzwiller value and the VMC estimate for weak and intermediate on-site 
correlation strengths U. We observe that the GWF gives a good estimate of 
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the correlation functions and that the errors in the VMC are not much larger. 
However, it is worth noting that the Gutzwiller method slightly underestimates 
the role of correlation as is seen from the higher than exact values of < niTij > 
and the lower than exact values of for | < SiSj > \ for both U = 1.0 and 
U — 4.0. The reduction in charge fluctauations and concommitant increase 
in spin fluctuations with increase in U is suggestive of the development of a 
spin at each site in the large U limit. In table (3), we present the exact, true 
Gutzwiller and VMC bond-orders for two different values of \J . We see that 
the VMC and the Gutzwiller bond-orders compare well with the exact values 
for wide range of correlation strengths. The larger bond order in the Gutzwiller 
method compared to the exact value is consistent with the observation that 
this approximation slightly underestimates the role of U . 

The results presented in tables (2) and (3) are obtained after averaging 
over equivalent site pairs and up and down spins. The VMC estimates before 
averaging are not strictly equal for equivalent pairs while the true Gutzwiller 
numbers reflect the symmetry of the system. 

The comparison between MO-VB and exact values in the ground state 
are favourable for threshold ~ bandwidth, although the MO-VB technique is 
designed to reproduce excitation gaps for various correlation strengths and is 
not the method of choice for the ground state properties. 
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3 Results on Fragments of Cqq 

In this section, we present our results on the fragments of Ceo using the MO- 
VE technique and compare these with results obtained from the VMC method. 
We have studied some catacondensed fragments of Cqq which contain charac- 
teristic chemical units of Cqq. These fragments are pyracylene, fluoranthene 
and corranulene, shown in fig(l). These fragments contain one or more five- 
memebered rings essential for displaying the minimal topological features of 
Ceo- The trends observed in these systems are expected to be indicative of the 
behaviour of the full molecule. 

3.1 Pyracylene 

The pyracylene molecule is a fragment of Ceo that is amenable to exact calcu- 
lations. The study of this molecule allows a further test on the approximations 
besides being of importance in its own right as the molecule is a characteristic 
fragment of Ceo. In table (4) we present the energy of pyracylene for vari- 
ous values of U. The energy obtained from the VMC calculation compares 
favourably with that from the exact calculation for U =1.0. 

In table (5a) we present the bond order data for the unique bonds in pyra- 
cylene. From the bond order data, we see that the 1-2 bond or its equivalents 
is the strongest bond in the system. This trend is found for all correlation 
strengths studied and is also seen in the exact PPP|^4| ground state. The 



external carbon atoms of the pentagons form almost independent ethylenic 
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units with a slight extended conjugation as seen from the rather small 2-3 
bond order. This is summarized in fig. 2 for U/t = 1.0. The bond orders of 
the internal bonds of the six membered rings are smaller than that of benzene 
for the same value of U (table(3)) although we can still view the electrons as 
being delocalized over the benzene rings. This picture is reinforced by the ring 
currents obtained in the PPP model which show large diamagnetic circulations 
around the naphthalenic perimeter but paramagnetic circulations around the 
pentagons. The similarity in PPP and Hubbard bond orders implies simi- 
lar ring currents in the Hubbard model also. We also expect paramagnetic 



circulations along the perimeter as found in the PPP model . 

The kinetic energy contribution to the total ground state energy decreases 
with increasing U and comparison with the PPP results allows us to estimate 
that the Hubbard model with an Ueff of nearly S.Ot should be similar to the 
PPP model. In table(5b) we present the bond order data of table(5a) normal- 
ized with respect to the 1-2 bond. Clearly, the Hubbard model shows smaller 
bond order variations compared to the PPP model. If the off-site interaction 
part in the PPP model is treated in the mean field, it would renormalize the 
transfer part of the Hamiltonian. It is for this reason that the interaction 
term in the PPP model infiuences the geometry, while the Hubbard interac- 
tions, which are purely on-site have a smaller influence on the geometry of the 
system. 

In tables(6a) and (6b) we present the spin and charge correlations for both 
on-site and nearest neighbour atom pairs and summarize them in flg. 2 for 
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U/t = 1.0. The 1-2 bond has a very large antiferromagnetic spin correlation 
(-0.402) although the sites 1 and 2 have rather large ionic contributions. In 
the half-filled case, the diagonal (on-site) spin and charge correlation functions 
in the large U limit are both 1.0, while in the f/ ^ limit, they are 0.5 and 
1.5 respectively.0 The large value oi < h\> and < sfsg > for the sites 1 and 
2 can be reconciled only if we assume that the spins are aligned antiparallel 
with a high probability whenever these sites are singly occupied. Thus the 
ethylenic units show no frustration for delocalization, a trend also evident 
from the bond order. The spin correlations follow the same trend as the bond 
orders, with spins at sites 2 and 3 being almost uncorrelated. This implies 
that the 2-3 bond is not frustrated. The most ionic sites in the molecule 
are the interior sites, 4 and 11. The exterior benzene sites are least ionic . 
The hexagon-hexagon bond is also found to be marginally stronger than the 
hexagon-pentagon bond, a feature in common with Cgo- 

3.2 Fluoranthene 

The dependance of the energy of fiuoranthene on U is given in table(4). Here, 
exact results are not available for comparison but the energy scales as expected 
with U . From the bond order data of fiuoranthene (fig. 3a), we see that the 

pentagon bonds are weak, while the bonds on the perimeter are strong. The 

^These results are obtained by computing 4(sfsf) and (nf), given that for the half-filled 
case, in the U/t = hmit, the probabihties of the Fock space states (|0), | t), | i)cmd\ tJ.)), 
at a given site i, are all equal to 0.25. In the large U limit, these probabilties are respectively 
0, 0.5, 0.5 and 0. 
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weakest bond is the bridging 4-5 bond. The bonds forming the benzenic unit 
and the naphthalenic unit are fairly strong. From the bond order data, we can 
visuahze fiuoranthene as weakly bridged benzene and naphthalene units, a fea- 
ture consistent with an earlier PPP model calculation. Furthermore, the bond 
order in the isolated hexagon is nearly the same as in benzene. The bonds 
in the naphthalenic unit show more alternation. We notice that even in fiu- 
oranthene, the hexagon-hexagon bond is slightly stronger than the hexagon- 
pentagon bond. The dependence of bond order on correlation strength is 
shown in table (7). All the bond orders decrease with increasinng U. The de- 
crease is sharper at larger U values. The bond orders tend to become uniform 
within the benzenic and naphthalenic units. 

The spin and charge correlations of the nearest neighbour bonds for U = 1.0 
are shown in fig(3b). The bridging 4-5 bond has a very small antiferromagnetic 
spin correlation, despite being a nearest neighbour bond. Thus we expect most 
of the frustration of the pentagon to be manifest in this bond. The charge- 
charge correlation is largest for this atom pair, and so is the on-site charge 
correlation. The small spin-spin correlation and the large charge correlations 
and site ionicities of this bond can be reconciled if we assume that both the sites 
are simultaneously empty or doubly occupied, and if singly occupied the spins 
are uncorrelated. It also imphes that the most reactive sites in fiuoranthene 
are likely to be these sites. In tables (8a) and (b) we present the dependence of 
charge and spin correlation son U. While all the sites tend to be more covalent, 
the relative ionicities as well as charge and spin correlations evolve smoothly 
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with U. We also observe that as U increases, the bond orders decrease, while 
the spin-spin correlations increase, a trend expected from the Hubbard model. 

3.3 Corannulene 

The energy of the ground state of corannulene for different values of U is pre- 
sented in table(4). The ground state energy per carbon atom, ec is the lowest 
for this fragment, with Sc for the fragments obeying the trend corannulene i 
fluoranthene j pyracylene. The stabilization decreases with increasing U, al- 
though the same trend persists. In fig. (4), we present the the bond order 
data of corannulene for U — 1.0. Unlike in the other two fragments of Cqq, the 
pentagonal unit in corannulene is completely enclosed. Therefore the bonds of 
the pentagon in this fragment most closely approximate the h-p bonds in Cqq. 
All h-p bonds in corannulene have almost equal bond orders. We observe that 
the h-p bonds are slightly weaker than the h-h bonds. The outer bonds of 
corannulene also fall into two categories. The 1-2 like bonds are stronger than 
the 2-3 like bonds. The higher bond orders of the former group suggests that 
bonds involving even a single three coordinated carbon atom tend to be weaker 
due to frustrations encountered during delocalization. The outer bonds with 
large bond orders evolve into the h-h bonds in Cgo- The bond orders show 
similar trends for different values of U (table (9)), although the magnitudes 
decrease with increasing correlation strength. 

The spin-spin and charge-charge correlations for [/ = 1.0 are shown in fig. 
(5). The interior sites are more ionic than the exterior sites. The charge-charge 
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correlations for atom pairs in the interior sites is larger than for the atom pairs 
in the exterior sites. The charge-charge correlation for the nearest neighbour 
pair involving one interior and one exterior atom is in between the values for 
the purely interior and the purely exterior atom pairs. The smallest charge- 
charge correlation is found for bonds with the largest bond order. The diagonal 
spin correlations are nearly uniform at all sites. The spin-spin correlations are 
nearly the same for any atom pair involving atleast one three connected carbon 
atom, while, for the atom pair with both two-connected carbons, it is much 
larger. This reflects the bond frustration in bonds involving three connected 
carbon atoms. These trends continue to hold for larger values of U, as seen 
from tables (10a) and (10b). 

4 VMC Ground State Properties of Cqq 

We have obtained the ground-state energy of neutral Ceo, with and without 
bond-alternation for several values of U. We have also computed the nearest 
neighbour bond orders and the charge and spin correlations of Cqo in both 
cases. The energies of Ceo are presented in table (11). These energy values are 



in good agreement with those of Krivnov et al for the values of U studied [p5 
We note that the energy per carbon atom is lower than that of any of the frag- 
ments studied. The electronic energy is further reduced on introducing bond 



alternation. However, we need a reliable estimate of the strain energy p6|, p7 | 



to be able to arrive at an equilibrium distortion of the bonds in Ceo. Nev- 



21 



ertheless, it is now established that the h-h bonds in Ceo are shorter than 



the h-p bonds[^, It is interesting to compare the energy per bond of 
the fragments with that of Cqq. The U = 1.0 values for benzene, pyracylene, 
fluoranthene, corranulene and Cgo are -1.10, -0.95, -0.99, -0.96 and -0.88 re- 
spectively. This trend is similar at higher values of U as well. This brings out 
clearly that Cqo is a more frustrated system than its fragments. The benzene 
value is clearly outside the range of the fragments of Ceo- 

In table (12) we present bond orders, spin-spin and charge-charge corre- 
lations for Ceo with and without bond alternation, for the h-h and h-p bonds 
as well as diagonal charge and spin correlations. We notice that the h-h bond 
orders are larger than the h-p bond orders for undistorted Ceo (uniform bond 
length) and the difference increases with alternation. While it is possible to 
optimise the bond length differences for the h-h and the h-p bonds on the ba- 
sis of Coulson's empirical formulae, we have not attempted it. From the bond 
orders for the uniform transfer integral, Coulson's formula]^, ^] predicts a 
h-h bond length of 1.429 A and a h-p bond length of 1.453 A. However, we 
have chosen the transfer integral of the h-h bond to be 20% stronger than that 
of the uniform case and the h-p bond to be 10% smaller than in the uniform 
case. This choice conserves the total transfer integral for the buckyball and 
allows comparison with the uniform case. The larger alternation has been 
employed to offset the errors in monte carlo estimates. 

The difference in bond order between the h-h and h-p bonds in uniform 
Ceo is larger than that in corannulene. This seems to indicate that the bond 
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frustrations in Ceo are more localized on the pentagons than in the fragments 
of Ceo or indeed any open frustrated systems. The h-p and h-h charge-charge 
correlations show the opposite trend when compared to corannulene. This 
is because the charge-charge correlations are smaller for large bond orders 
as the transfer of an electron requires the dissimilar occupancy of the sites 
between which hopping is occuring. The charge-charge correlation for the h- 
h bond is smaller than for the h-p bond in Ceo- The spin-spin correlation 
is more antiferromagnetic for the h-h bond than for the h-p bond in Ceo- 
The h-h bond in Ceo has weaker antiferromagnetic fluctuations than the h-h 
bond in corannulene, while the h-p bond in Ceo has smaller antiferromagnetic 
fluctuations than the h-p bond in corannulene. This trend holds for all values 
of U. 

The site ionicities of the carbon sites in Ceo are smaller than the interior 
sites in corannulene, as can be seen from the smaller diagonal charge correlation 
in the former. This is also confirmed by the larger diagonal spin correlation in 
Ceo compared with that for the pentagonal sites in corannulene. 

The spin and charge correlations for the h-h bond increase when alternation 
in the transfer integral is introduced, while these quantities decrease for the 
h-p bond. The site ionicities are smaller than those of uniform Ceo. This 
trend persists for all values of U studied by us. 

The average bond order of a molecule gives the extent of kinetic stabiliza- 
tion due to delocalization. It is interesting to compare the average bond orders 
for U/t — for benzene. Ceo and graphite[?], which are 2/3, 0.518 and 0.525. 
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In benzene and Ceo the bond orders are marginally reduced for U/t — 1.0 
to 0.66 and 0.513 respectively. A similar trend is expected in graphite, upon 
introducing correlations. Thus, it appears that in Cqq the kinetic stabilization 
is only slightly smaller than in graphite, in spite of bond frustration in the 
former due to the five-membered rings. This result could be an artefact of 
the Hubbard model; for even in the mean-field limit, extended range interac- 
tions renormalize the transfer integrals, leading to contributions to the kinetic 
energy and hence to the bond orders. It is also interesting to compare the pen- 
tagons in Ceo with the cyclopentadienyl radical. The pentagonal bond order 
in this molecule is 0.58 while the h-p bond order in Ceo is 0.47 for U — 1.0. 
Thus the pentagons in Ceo are more frustrated than a single pentagonal unit. 

The value of the Gutzwiller parameter g for which the energy is a minimum 
indicates the overall extent of supression of the doubly occupied sites in a given 
system. In table (13), the value of the Gutzwiller parameter g for which the 
energy is a minimum is given for all the systems reported in this paper. There 
is no systematic variation in g across the fragments and all the systems seem 
to suppress double occupancies to a similar extent for a given value of U . 

5 Summary 

We have studied the fragments of Ceo by various CI techniques as well as 
the VMC. Comparison of these studies allows us to derive reasonable monte 
carlo parameters for the simulation of larger systems hke Ceo. The bond 
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order analysis of these molecules allows us to visualise structures which are 
consistent with the spin and charge correlations. Wc are able to observe the 
evolution of the properties of the fragments towards the properties of Ceo 
through the properties of the hexagon-hexagon like and hexagon-pentagon like 
bonds. Corannulene is found to be the fragment of Ceo that approximates the 
full molecule most closely. We have studied the uniform and bond alternated 
models for Ceo- Even when the transfer integral is uniform, the bond orders 
of Ceo show more marked alternation than the fragments. 

Even though the VMC is an approximation over and above that inherent 
in the Gutzwiller wavefunction it is interesting to see that the VMC technique 
provides a reasonable description of the ground state electronic structure of 
molecular systems. The VMC energies are in good agreement with model 
exact energies for intermediate values of U/t, but progressively deteriorates 
for large U/t. We are currently involved in extending the Gutzwiller type of 
VMC approach to models with nearest neighbour Coulomb interactions. 
Acknowledgements: We thank Y. Anusooya for assistance rendered in the 
MO-VB computations on PPP models and Dr. Biswadeb Dutta for help with 
the computational facility at JNCASR. 
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Table 1. Ground state energy of benzene (in units of as a function of U 

from exact and approximate calculations. The Qmin values are given in 
parantheses. The VMC energies do not change on doubling the number of 





MPE for sampling 


u 


Eexact 


Emo-vb 


Eoutz 


EvMC 


0.5 


-7.2744 


-7.2742 


-7.2741 (0.91) 


-7.274 (0.91) 


1.0 


-6.6012 


-6.5963 


-6.5964 (0.83) 


-6.597 (0.83) 


2.0 


-5.4105 


-5.3744 


-5.3868 (0.69) 


-5.388 (0.71) 


4.0 


-3.6697 


-3.3516 


-3.5371 (0.47) 


-3.538 (0.46) 



Table 2. Spin-Spin correlations (4 < s|s| >) and charge-charge correlations 

(< niUj >) of benzene for U — 1.0 and 4.0, from exact and approximate 
calculations. The VMC values do not change on doubling the Monte Carlo 

sampling size. 



u 


hJ 


4 < 


> 


< niUj > 






Exact 


Gutz 


VMC 


Exact 


Gutz 


VMC 


1.0 


1,1 


0.567 


0.565 


0.557 


1.432 


1.435 


1.443 


1.0 


1,2 


-0.267 


-0.265 


-0.245 


0.816 


0.817 


0.812 


1.0 


1,3 


0.022 


0.016 


0.011 


0.989 


0.989 


0.976 


1.0 


1,4 


-0.077 


-0.068 


-0.062 


0.960 


0.953 


0.952 


4.0 


1,1 


0.778 


0.751 


0.750 


1.222 


1.249 


1.250 


4.0 


1,2 


-0.435 


-0.404 


-0.395 


0.905 


0.911 


0.911 


4.0 


1,3 


0.140 


0.095 


0.077 


0.999 


0.979 


0.979 


4.0 


1,4 


-0.188 


-0.131 


-0.113 


0.989 


0.971 


0.971 
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Table 3. Bond-orders of the nearest neighbour bond in benzene for U = 1.0 
and 4.0 from exact and approximate calculations. The VMC values do not 
change on doubling the Monte Carlo samphng size. 







u 


Exact 


Gutz 


VMC 


0.5 


0.6646 


0.6646 


0.664 


1.0 


0.6582 


0.6584 


0.660 


2.0 


0.6316 


0.6344 


0.636 


4.0 


0.5278 


0.5441 


0.545 



Table 4. Ground state energies of pyracylene, fluoranthene and corannulene 
(in units of t) from restricted CI and VMC calculations. Numbers in 
parantheses are the energies per site. MO-VB energies are not given for 
larger values of U as the technique is known to give poor ground state 
energies for large U. The VMC values do not change on doubling the Monte 

Carlo sampling size. 







u 


Pyrac-ylono 


Fluorantheue 


Coranuulene 




Emo-vb Evmc 


Emo-vb Evmc 


Emo-vb Evmc 


0.5 
1.0 
2.0 
4.0 


-17.7002 -17.713(-1.27) 
-16.0631* -16.116(-1.15) 
-12.9817 -13.235(-0.95) 
— -8.666 (-0.62) 


-20.5231 -20.558(-1.28) 
-18.5943 -18.730(-1.17) 
-14.8631 -15.422(-0.96) 
— -10.164 (-0.64) 


-26.2689 -26.306(-1.32) 
-23.8655 -24.007(-1.20) 
-19.2221 -19.847(-0.99) 
— -13.158(-0.66) 



* exact energy obtained from full CI calculation in the singlet subspace of 
dimension 2760615 is -16.13 
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Table 5a. Bond-orders of the unique bonds and total kinetic energy (KE) in 
the ground state of pyracylcnc from VMC calculations. The VMC values do 
not change on doubling the Monte Carlo sampling size. 





C/ = 0.5 


U = 1.0 


C/ = 2.0 


C/ = 4.0 


PPP (exact) 


1,2 


0.819 


0.819 


0.795 


0.695 


0.816 


2,3 


0.399 


0.402 


0.382 


0.330 


0.239 


3,4 


0.517 


0.509 


0.494 


0.442 


0.455 


5,6 


0.652 


0.647 


0.630 


0.562 


0.697 


6.9 


0.609 


0.605 


0.581 


0.518 


0.604 


4,11 


0.531 


0.531 


0.519 


0.163 


0.526 


KE 


-19.365 


-19.228 


-18.594 


-16.448 


-17.860 
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Table 5b. Normalized bond orders of the unique bonds in the ground state of 



pyracylene from VMC calculations. 







hj 


t/ = 0.5 


U = 1.0 


C/ = 2.0 


C/ = 4.0 


PPP (exact) 


1,2 


1.000 


1.000 


1.000 


1.000 


1.000 


2,3 


0.490 


0.491 


0.480 


0.474 


0.293 


3,4 


0.631 


0.621 


0.621 


0.636 


0.558 


5,6 


0.795 


0.790 


0.793 


0.809 


0.854 


6,9 


0.790 


0.739 


0.731 


0.745 


0.740 


4,11 


0.618 


0.618 


0.652 


0.665 


0.615 
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Table 6a. On-site and nearest neighbour VMC spin-spin correlations 
(4 < s|sj >) of pyracylene for different U /t values. The VMC values do not 
change on doubling the Monte Carlo sampling size. Numbers in parantheses 
are the exact AO-VB results, given for comparison. 





Nearest nei{ 


^hbour 






On- 


-site 






U = 1.0 


U = 2.0 


U = 4.0 


hj 


U = 1.0 


U = 2.0 


U = 4.0 


1,2 


-0.402 (-0.4308) 


-0.475 


-0.595 


1,1 


0.565(0.5849) 


0.625 


0.744 


2,3 


-0.099 (-0.0930) 


-0.107 


-0.138 


3,3 


0.551(0.5563) 


0.616 


0.711 


3,4 


-0.155 (-0.1532) 


-0.182 


-0.241 


4,4 


0.549 (0.5477) 


0.604 


0.706 


5,6 


-0.250 (-0.2598) 


-0.254 


-0.382 


6,6 


0.554 (0.5651) 


0.617 


0.719 


4,11 


-0.164 (-0.1703) 


-0.205 


-0.273 










6,9 


-0.225 (-0.2219) 


-0.260 


-0.326 
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Table 6b. On-site and nearest neighbour VMC charge-charge correlations 
(< niTij >) of pyracylene for different U/t values. The VMC values do not 
change on doubling the Monte Carlo sampling size. 





Nearest 


neighbour 




On-site 


I- J 


U = 1.0 


U = 2.0 


U = 1.0 


'■J 


U = l.C 


U = 2.0 


U = 1.0 


1,2 


0.794 


0.841 


0.929 


1,1 


1.514 


1.444 


1.321 


2,3 


0.964 


0.972 


0.992 


3,3 


1.431 


1.374 


1.281 


3,4 


0.964 


0.984 


1.012 


4,4 


1.621 


1.570 


1.400 


5,6 


0.737 


0.776 


0.846 


6,6 


1.228 


1.240 


1.171 


6,9 


0.702 


0.738 


0.810 










4,11 


1.059 


1.064 


1.393 
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Table 7. Bond-orders of the unique bonds and total kinetic energy (KE) in 
the ground state of fluoranthene from VMC calculations. The VMC values 
do not change on doubling the Monte Carlo sampling size. 





C/ = 1.0 


C/ = 2.0 


C/ = 4.0 


1,2 


0.643 


0.612 


0.518 


2,3 


0.666 


0.637 


0.550 


3,4 


0.611 


0.591 


0.519 


4,5 


0.388 


0.377 


0.320 


5,6 


0.657 


0.630 


0.548 


6,7 


0.611 


0.593 


0.514 


7,8 


0.707 


0.676 


0.579 


8,9 


0.546 


0.522 


0.451 


4,14 


0.566 


0.539 


0.472 


5,16 


0.512 


0.494 


0.436 


9,16 


0.526 


0.512 


0.446 


KE 


-22.256 


-21.397 


-18.539 
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Table 8a. On-site and nearest neighbour VMC spin-spin correlations 
(4 < sfSj >) of fluoranthene for different U/t values. The VMC values do not 
change on doubling the Monte Carlo sampling size. 



Nearest neighbour 


On-site 


hJ 


U = 1.0 


U = 2.0 


U = 1.0 


i-J 


U = 1.0 


U = 2.0 


U = 4.0 


1,2 


-0.259 


-0.300 


-0.395 


1,1 


0.5617 


0.6358 


0.7548 


2,3 


-0.284 


-0.327 


-0.409 


3,3 


0.5702 


0.6372 


0.7399 


3,4 


-0.228 


-0.283 


-0.362 


4,4 


0.552 


0.619 


0.729 


4,5 


-0.087 


-0.109 


-0.137 


5,5 


0.561 


0.622 


0.729 


5,6 


-0.260 


-0.305 


-0.403 


6,6 


0.552 


0.621 


0.738 


6,7 


-0.222 


-0.271 


-0.360 


7,7 


0.565 


0.630 


0.745 


7,8 


-0.303 


-0.361 


-0.452 


8,8 


0.557 


0.626 


0.744 


8,9 


-0.176 


-0.206 


-0.280 


9,9 


0.551 


0.620 


0.731 


4,14 


-0.187 


-0.220 


-0.298 


16,16 


0.561 


0.624 


0.722 


5,16 


-0.167 


-0.211 


-0.271 










9,16 


-0.166 


-0.176 


-0.233 
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Table 8b. On-site and nearest neighbour VMC charge-charge correlations 
(< riirij >) of fluoranthene for different U/t values. 





Nearest 


neighbour 






On 


-site 






U = 1.0 


U = 2.0 


= 4.0 


ij 


U = 1.0 


U = 2.0 


f/ = 4.0 


1,2 


0.835 


0.874 


0.926 


1,1 


1.452 


1.373 


1.249 


2,3 


0.807 


0.859 


0.911 


3,3 


1.397 


1.364 


1.260 


3,4 


0.864 


0.898 


0.945 


4,4 


1.523 


1.433 


1.318 


4,5 


1.013 


1.003 


1.016 


5,5 


1.515 


1.443 


1.321 


5,6 


0.800 


0.841 


0.909 


6,6 


1.344 


1.291 


1.209 


6,7 


0.790 


0.825 


0.883 


7,7 


1.435 


1.373 


1.234 


7,8 


0.839 


0.849 


0.861 


8,8 


1.379 


1.294 


1.193 


8,9 


0.839 


0.860 


0.914 


9,9 


1.452 


1.384 


1.278 


4,14 


0.934 


0.939 


0.981 


16,16 


1.508 


1.452 


1.343 


5,16 


0.965 


0.986 


0.998 










9,16 


0.920 


0.949 


0.981 
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Table 9. Bond-orders of the unique bonds and total kinetic energy (KE) in 
the ground state of corannulene from VMC calculations. 





[/ = 1.0 


C/ = 2.0 


C/ = 4.0 


1,2 


0.745 


0.721 


0.641 


2,3 


0.529 


0.509 


0.454 


16,17 


0.504 


0.489 


0.433 


3,17 


0.551 


0.535 


0.491 


KE 


-28.562 


-27.673 


-24.751 



Table 10a. On-site and nearest neighbour VMC spin-spin correlations 
(4 < s^Sj >) of corannulene for different U/t values. 



Nearest neighbour 


On-site 


hJ 


[/ = 1.0 


C/ = 2.0 


C/ = 4.0 


hJ 


U^l.O C/ = 2.0 


C/ = 4.0 


1,2 


-0.323 


-0.380 


-0.411 


1,1 


0.552 0.612 


0.720 


2,3 


-0.155 


-0.194 


-0.240 


3,3 


0.539 0.604 


0.702 


16,17 


-0.148 


-0.183 


-0.212 


16,16 


0.537 0.601 


0.699 


3,17 


-0.173 


-0.211 


-0.278 
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Table 10b. On-site and nearest neighbour VMC charge-charge correlations 
(< riiTij >) of corannulene for different U/t values. 





Nearest 


neighbour 






On 


-site 






U = 1.0 


U = 2.0 


U = 4.0 


^■J 


U = 1.0 


U = 2.0 


U = 4.0 


1,2 


0.725 


0.777 


0.847 


1,1 


1.414 


1.359 


1.251 


2,3 


0.840 


0.870 


0.915 


3,3 


1.422 


1.357 


1.284 


16,17 


0.997 


1.001 


1.017 


16,16 


1.567 


1.494 


1.378 


3,17 


0.895 


0.917 


0.957 











Table 11. Ground state energy of uniform and bond alternant Ceo for various 
values of U from VMC. 



u 


uniform 


bond alternant 




E 


E/site 


E/bond 


E 


E/site 


E / bond 


1.0 


-78.891 


-1.315 


-0.877 


-81.448 


-1.358 


-0.905 


2.0 


-66.074 


-1.101 


-0.734 


-68.733 


-1.146 


-0.764 


3.0 


-55.669 


-0.928 


-0.619 


-57.473 


-0.958 


-0.639 


4.0 


-44.678 


-0.745 


-0.496 


-47.768 


-0.796 


-0.531 


5.0 


-36.088 


-0.602 


-0.401 


-39.566 


-0.695 


-0.440 


6.0 


-28.787 


-0.480 


-0.320 


-32.435 


-0.541 


-0.360 
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Table 12. On-site and nearest neighbour VMC charge-charge correlations 
(< niUj >) and spin-spin correlations (4 < >) and bond orders of 
uniform and bond alternant Ceo for different U/t values. 



u 




uniform 


bond alternant 






h-h 


h-p 


diag 


h-h 


h-p 


diag 


1.0 


bo 


0.597 


0.471 




0.849 


0.367 






s-s 


-0.212 


-0.132 


0.550 


-0.296 


-0.099 


0.5500 




c-c 


0.817 


0.901 


1.150 


0.781 


0.927 


1.150 


2.0 


bo 


0.582 


0.461 




0.827 


0.359 






s-s 


-0.241 


-0.152 


0.589 


-0.336 


-0.115 


0.601 




c-c 


0.872 


0.919 


1.402 


0.819 


0.938 


1.399 


3.0 


bo 


0.535 


0.414 




0.756 


0.321 






s-s 


-0.314 


-0.191 


0.690 


-0.431 


-0.143 


0.699 




c-c 


0.912 


0.947 


1.310 


0.879 


0.960 


1.301 


4.0 


bo 


0.448 


0.336 




0.634 


0.254 






s-s 


-0.404 


-0.233 


0.787 


-0.545 


-0.169 


0.799 




c-c 


0.949 


0.970 


1.258 


0.928 


0.977 


1.201 
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Table 13. Minimum value of the Gutzwiller parameter g for fragments of Cqq, 
uniform and bond alternant Ceo (A and B). 



u 


in'rac-yleiie 


fluoraiitlioue 


corauiiuloiio 


aK,(A) 


Coo(B) 


0.5 


0.91 


0.92 


0.93 






1.0 


0.85 


0.83 


0.88 


0.85 


0.87 


2.0 


0.69 


0.69 


0.71 


0.73 


0.74 


4.0 


0.50 


0.47 


0.51 


0.53 


0.53 
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Figure Captions 



Figure 1. Structures of Pyracylene, Fluoranthene and Corannulene. 

Figure 2. Bond orders(A), diagonal charge-charge (spin-spin) correlations(B), 
narest neighbour charge-charge correlations(C) and nearest neighbour 
spin-spin correlations (D) for unique sites and bonds of pyracylene for 
U = 1.0. 

Figure 3a. Bond orders and diagonal charge-charge (spin-spin) correlations 
for unique bonds and sites in fluoranthene for U — 1.0. The axes of 
symmetry are shown. 

Figure 3b. Unique nearest neighbour charge-charge correlations and spin- 
spin correlations of fluoranthene for U — 1.0. 

Figure 4. Bond orders and diagonal charge-charge (spin-spin) correlations 
of in corannulene for U — 1.0. The difference between symmetry related 
quantities indicates the extent of scatter in VMC estimates. 

Figure 5. Nearest neighbour charge-charge correlations and spin-spin cor- 
relations in corannulene for U = 1.0. The difference between symmetry 
related quantities indicates the extent of scatter in VMC estimates. 
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